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This presentation briefly describes the Bayesian Block algorithm in the context of its apphcation to analysis of 
time series data from the Fermi Gamma Ray Space Telescope. Details are to be found in a recently published 
paper. 
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1. Introduction 

Much astrophysical information obtained with the 
Fermi Gamma Ray Space Telescope arises from analy- 
sis of time series data. Accordingly some of the desired 
capabilities of science analysis software are the ability 
to: 

• Detect transient events, such as Gamma Ray 
Bursts, and flares in active galactic nuclei and 
other variable objects such as the Crab Nebula. 

• Characterize these events by modeling the shape 
of the light curve. 

• Estimate values of parameters from this char- 
acterization, such as rise and decay times, time 
scales of variability, and an index of variability. 

• Identify times over which the near constancy of 
the flux suggests good time intervals over which 
to carry out estimates of the energy spectrum 



2. The Bayesian Blocl( Algorithm 

The Bayesian Blocks algorithm provides a useful 
method to address these problems. The concept un- 
derlying this approach is to find the best possible 
representation of time- sequential data as a series of 
segments over which the underlying signal is con- 
sta nt to within the o bservational errors. The pa- 
per [Scargle et al.ll2013| describes a dynamic program- 
ming algorithm which yields such an optimal price- 
wise constant model for time series data in general, 
and time-tagged photon data as a special case. 

All operations of Bayesian Blocks can be imple- 
mented with no limitation on sampling, time resolu- 
tion or signal amplitude. Data gaps and variable ex- 
posure are easily accommodated. Other applications 
include data-adaptive histograms, multi-variate time 
series analysis and a novel approach to delay estima- 
tion in strong lensing events. 



3. Published Code 



De tails and principles are given in [Scargle et al.l 
|2013| . This reference includes data files and Matlab® 
code for implementing the Bayesian Blocks algorithm 
and reproducing the figures in the paper, thus im- 
plementing the discipline of Rep roducible Research 
as invented in IClaerboutI Il990j | and devleoped in 
[Donohoet al.ll2008l| . 

Here are a few corrections to this paper and the 
attached code: 

• In the text of the A strophysical Journal paper 
Scargle et al.ll2013{ the natural logarithm sym- 



bol was omitted from Equation 21, which should 
be 

ncp_prior = 4 - log(73.53poA^"'^'^*) 

This error was reported by Peter Williams, at 
the Center for Astrophysics, Harvard University. 
The attached code implements this expression 
correctly. 

On hue 222 of MatLab script 
named f ind_blocks.m the expression 
nuin_this = sum( nn_vec( ii_2: ii_2 ) ); 
should be 
nuin_this = sumC nn_vec( ii_l: ii_2 ) ); 



That is, the summation of the populations of all 
of the data cells in the block is over the index 
range from ii_l to ii_2. This error, only af- 
fecting the display of the block model and not 
its computation, was reported by Mike Wheat- 
land, the University of Sydney. 

• In the same code package, the first line of the 
function script plot_blocks_meas .m had a 
typo in the function name. It should read 
function rate_vec = plot_blocks_meas( 
cpt_times, tt_vec, xx_vec, er_vec ). 
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In most implementations of MatLab this error 
causes no problems. 

Also it should be noted that this A-IatLab code was 
developed with an earlier version (7.9, from 2009) and 
there may be minor compatibility issues with more 
recent versions. As additional updates or corrections 
to this material are developed they will be posted at 
http : //bayesianblocks . blogspot . com 

In addition Jake Vanderplas'related blog "Dynamic 
Programming in Python: Bayesian Blocks" nicely 
described the algorithmic approach of dynamic 
programming, with examples in the context of 
histograms: http : // j akevdp . github . io/blog/ 

2012/09/ 12/dyncmiic-progrELmming-in-pythoii/. 
See also the posting at the Starship Asterisk dis- 
cussion forum: http://asterisk.apod.com/ 
viewtopic . php?f =35&t=29458 
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